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slowly  with  position  and  time 
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UyV  .VKI1  me, 'hods  we  have  considered  'Hie  propagation  of  electromagnetic  waves  In 
isotropic  lossless  media  which  vary’  slowly  wl(t>v  both  position  and  time.  It.  is  found  that  in 
such  media  the  meaning  of  various  quantities,  such  tis  the  group  velocity,  must  be  reinterpreted. 
Tho  theory  is  applied, \to  study  propagation  in  space-time,  varying  dielectrics,  and  plasmas. 


I.  INTRODUCTION 

A-diflicult  problemtoconsidcrcoaceptunllyiswavc 
propagation  in  media  which  vary  with  botlt  position 
and  time.  In  such  media  our  standard  concepts  of 
frequency,  wavenumber,  and  group  velocity  no  longer 
apply.  Thai  is,  wc  find  that  frequency  can  be  defined 
only  as  the  tints  derivative  of  the  phase  function  and 
that  the  quantity  so  defined  is  a  function  of  both  the 
position  and  the  time  at  which  the  observation  is 
made.  Wc  also  find  that  the  group  velocity  no  longer 
retains  its  conventional  meaning,  in  spatially-homo- 
gcncous,  time-invariant  media  the  group  velocity  is 
interpreted  as  the  velocity  at  which  wave  packets 
centered  around  some  wavevcctor  k,  propagate 
[Jeffreys  and  Jeffreys ,  1962].  In  space-time  varying 
media  this  is  no  longer  true.  In  fact,  as  we  shall  sec, 
values  of  and  k  do  not  propagate  with  the  group 
velocity  V»w;  it  is  rather  different  quantities  (which 
arc  functions  of  w  and  k)  which  arc  propagated  at 
this  velocity.  Tire  same  conclusions  hold  true  for 
energy  (low  (i.c.,  the  energy  flux  does  not  propagate 
with  the  group  velocity). 

In  this  paper  we  will  study  the  propagation  of 
electromagnetic  waves  in  lossless  media  which  vary 
slowly  with  (position  and  time.  Wc  will  therefore 
employ  the  four-dimensional  WKB  (Wentzcl* 
Kramcrs-Brillouin)  method,  The  WKB  method  was 
first  applied  in  three  dimensions  by  Sommerfeld  and 
Range  [191 1],  That  is,  in  media  in  which  the  properties 
depended  oil  position  x  but  not  on  time,  Sommerfeld 
and  Rungc  considered  solutions  of  the  form 
exp  (Lt  j"  k-r/x).  A,  a  consequence  of  the  fact  that 
the  phase  function  must  be  uniquely  defined,  Sommer- 
fold  and  Runge  then  concluded  that  V  x  k  »  0  was 


required.  This  is  the  original  version  of  thc  Sommcr- 
feW- Range  jaw.  This  result  can  be  cxtcndcdjto  the 
four-dimensional  case,  as  has  been  done  by  WhUliam 
[1960]  and  Poecerlein  [1962],  That  is,  in  space-time 
varying  media  the  existence  of  a  uniquely  defined 
wave  function,  exp  (/  !  («  dt  -  k*r/x)]  requires  that 
Vo  —  flk/j )/  <=>  0,  This  is  the  four-dimensional 
Sommerrcld-Rimgc  jaw. 

In  the  present  paper  wc  will  first  give  an  elementary 
derivation  of  tiie  four-dimensional  Sommcrfeld-Rungc 
law.  We.  will  then  examine  ils  implications,  and 
imaHy  indicate  its  use  in  studying  electromagnetic 
wave  propagation  in  isotropic,  lossless  media:  which 
vary  slowly  with  position: and  time. 

2.  GENERAL  THEORY 

2.1.  Discussion  of  the  Generalized  Sommerfeld- Range 
Law 

In  a  medium  which  varies  slowly  (c.g.,  in  a  dielectric 
the  conditions  for  slow  variation  arc  that  k3 »  V*k, 
U3  »  Oce/dt,  co  »  a/a/(lnt),  and  k  »  (l/<)|Vt|, 
where  e  is  the  permittivity  of  the  dielectric)  with  botii 
position  and  time,  the  WKB  approximation  for  the 
electric  field  strength  can  be  written  ns 

E(x,  I)  =>  e,(x,  0  exp  J  (k ulx  —  w  </r)J  (1) 

wlicre  e,(x,  t)  varys  slowly  witli  position  and  time 
compared  with  the  exponent,  and  L  is  *v  line  integral 
in  four  space  between  some  initial  point  (x„  lt)  and 
the  observation  point  (x,  t),  flic  functions  w  and  k 
are  generally  related  through  a  dispersion  relation 
of  the  form 
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w  «  ll'(k,  x,  I) 


(2) 
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or 

k  “  A'(u,  x,  i)  (3) 

where  k  =»  |k|.  For  example,  in  a  plasma  we  have 
k  =*•  —  «,J(x,r)]l/J,  where  c  is  the  speed  pf 

light,  and  w,  is  the  electron  plasma  frequency.  Now 
since  the  exponent  in  (I)  represents  a  phase  function, 
its  value  must  be  unique.  That  is,  the  line  integral  & 
should  be  independent  Of  the  path,  which  implies 
that  the  integrand  is  a- perfect  differential.  Therefore, 
there  must  exist  some  function  <£,  such  tluitk  = 
and  «  =•  —04>/0t.  As  a  consequence  of  these  relations 
we  then  have 

ilk/ dt  -|-  Vw  a  0  M) 

Equation  4  is  the  fo.'Mlimensionat,  version,  .of  the 
SommerfekURunge  law.  This  result  can  tuternutsly 
be  derived  by  employing  results  oij-wavo  cQnscrvakQtlj; 
as  has  been  done  by  Lighthlll  anil  IVhltham  [1055]. 
Equation  4  has  also  been  considered  by  WiUhnm 
[1960],  Pocvcrlcln  [1962],  and  Landau  and  Lfsehltz 
[1959]. 

2.2.  Discussion  of  thc<  Properties  of  Equation  4 

2.2.1.  Generalized  group  velocity.  To  discuss  the 
properties  of  (4)  in  a  medium  which  varies  with  both 
position  and  time,  let  us  substitute  (2)  into  (4),  Using 
the  fact  that  V  xk  »  0,  we  get 

0k/£M*(V*V)k»  “(VUOm  15) 

where 

The  quantity  V  In  (6)  can  be  interpreted  us  a  general* 
feed  group  velocity  as  we  shall  see  in  the  following 
discussion.  Upon  defining  d/dt  «  O/Ot  •}*  V> V 
equation  5  can  be  rewritten  as 

dk/dt  -  -(VfiOk.,  (7) 

From  (7)  we  see  that  if  }V  docs  not  depend  explicitly 
on  position,  then  dk/dt  »  0.  This  means  that  if  one 
moves  along  the  ray  with  the  velocity  V,  he  will 
observe  com'ant  values  of  k  (i.e.,  the  observer  moving 
With  V,  measures  constant  values  of  wavelength). 
Therefore,  if  IF  is  a  function  of  k  and  t,  but  -iocs  not 
depend  explicitly  on  x,  then  k  is  a  constant  of  the 
motion,  so  that  wave  packets  sharply  centered  mound 
some  wavenumber  k*  will  be  propagated  with  V  -• 

0  W/Ok  evaluated  at  k  *=  k„. 

To  examine  the  other  limit  when  tne  properties  of 


the  medium  vary  with  position  but  not  with  time, 
.it  is  convenient  ui  use  (3)  in  (4).  For  this  ease  we  obtain 
'(assuming  the  medium  is  isotropic) 

Ou/Ot  H-  V-Vic  -  ,  (8) 

From.  (8)  it  is  clear  that  if  A'  depends  on  w  and  x,  but 
does  not  explicitly  .depend  on  time,  then 

du/dt  “  0  (9) 

This,  means  that  the  observer  moving  along  the  ray 
with  the  velocity  V  will  observe  constant  values  of « 
(i.e.,  constant  wave  period).  Snell’s  law  follows 
imnicdiatcly  from  (9).  Therefore  wave  packets  sharply 
centered  about  some  frequency  wo  will  propagate  with 
the. group  velocity  V  (evaluated  at  «  «  w„). 

For  the  general  ease  when  the  properties  of  the 
medium  depend  on  hotly  .position  and  lime,  neither 
to  nor  k  wilt  be  invariant  ns  one  moves  along  the  ray 
with  the  group  velocity  V,  (However,  hi  the  next 
subsection  we  do  show  that  there  are  velocities 
Vlk>  V  and  V<w5  V  with  which  values  of  k  and  id 
nre  propagated  in  space-time  varying  media).  In  this 
ease  other  quantities  will  be  invariants  of  the  motion. 
Consider  each  scalar  component  of  (5),  NVc  have 

a  ii' 

Dk,/0t  +  (V*  V)A’(  -  -Jg-  GO) 

To  solve  (10)  we  consider  the  subsidiary  set 

dt  «*  f/.Vi/Vi  **  dxj  IT  dxj  IT  «  -  dk,/(0  »'/£>*,) 

{11V 

Let  ms  denote  the  pnrtienlar  iivtegrals  of  (I!)  by 
f(ku  NT  0  -  C„  g(kh  x,  t )  «  C„  bfku  x,  t)  -  C„ 
\f(k„  x,  I )  w  C,.  Then  it  can  be  shown  that  [Sneddoft, 
1957]  (a)  the  general  solution  of  (10)  Is  given  by 

C,  -  <I>(C„  Co  C,)  (12) 

where  4*  Is  an  arbitrary  function,  determined  by  the 
boundary  conditions  imposed  and  (b)  the  invariants 
of  motion,  for  an  observer  moving  with  the  velocity  V 
given  by  (0),  nre  C„  C\,  C„  and  CV  That  is,  each 
C,  satisfies 

dC,/dt  (.r,«M  V*V-»C.  O  113) 

Therefore,  for  medi »  which  vary  with  both  po  Uioft 
and  time,  the  (.na.ili/ed  group  vehvitv  Alined 
in  (6)  is  the  velocity  with  which  the  quantifies  i  ,  are 
propagated.  U  1'  only  indite  limit  »<f  tHy  homo* 
getieoiw,  limcdnvjiimt  uu.tn  tint  t.om,:  *4  the 
constant s-t’,  can  be  i.  Unified  with  an  .  k 


In  section  3.1  we  will  calculate  the  constants  of 
motion  for  several  examples  of  media  which  vary 
with  position  and  time. 

2.2.2.  The  angle  between  V  anil  k,  in  this  sub¬ 
section  we  will  demonstrate  that  in  isotropic  lossless 
media  the  normal,  k  -  (k/A*),  to  the  pluise  surface 
/)  lies  in  tiie  same  direction  as  the  group  velocity 
V.  NVc  will  also  calculate  the  angle  between  k  and  V 
for  the  case  of  a  simple  uuisotropy. 

fet  us  consider  (3)  for  the  case  wherftlic  dispersion 
region  can  be  written  as 

A-  «  A'(b,  0,  x,  s)  (14) 

v-hcre  0  |s  the  angle  which  the  vector  k  niukcs  with 
tlifc*  axis,  (Tbwcorrespcmds  to  the  case  in  which  the 
dispersion  relation  is  given  by  the  Appclton-llartrcc 
formula),  Using  (14)  in  '/>),  we  find,  upon  difleren- 
-tiating  impivltiy^tliat  the-group-velocityis- 

V  4r  <0w/i>A)[k  -  (6/k)(0  K/0  0)]  (15) 

Upon  tuking  the  del  and  cross  products  of  k  with  (1$) 
we  obtain  for  the  angle  y  between  k  and  V 

linn  y\  m  |k  «  V|/k«  V  •-  A'1  ~~  «  A“  |f  (Id) 

For  isotropic  media,  (k  independent  of  0)  the  angle 
7-0.  liquation  16  is  a  well-known  result  from  tiie 
Study  of  Whistlers  [Halt  and  Haskell,  1965;  Kelso, 
1964),  We  imve  shown  here  that  the  result  also  applies 
to  media  which  vary  slowly  with  both  position  and 
time,  provided  we  understand  that  V  Is  not  tiie  velocity 
with  which  values  of  k  and  «  arc  propagated,  One 
final  result  of  interest  Is  to  use  (16)  and  (15)  to  calculate 
the  magnitude  of  the  group  velocity.  We  get' 

ii-i  *»  w/'Umvr'  07) 

2.21  The  equation  of  motion  ofU.  Ill  this  subsection 
we  shall  demonstrate  that  the  classical  equation 
presented  by  landau  and  llfsehit:  [1959]  for  the 
motion  of  the  normal  to  <tfx,  /),  can  be  generalized 
to  include  space-time  varying  media.  For  isotropic 
media  we  write,  using  k  *  Akin  (5) 

AC dtt  dl)  I  Hulk  dll  -  -tVH'hi  08) 

Wo  r>cXt  wvo  the  di-perMou  uf  ilon  of  (1  in  the 
second  term  on  the  left-hand  side  of  US).  We  obtain, 
alter  r,  grouping  team 

»ufk  m  t  fa*  ,/m 

{  <•«*.. ,  1  V'iVA»  ,1  iV»fe  ,  1 19) 

where  ¥  V  <1  ,«„»  *  Mo*  substitute  (8)  i,»r 

(h  ntt  »UU)  T1-.  n  »!(  h 
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k(dk/dt)  -Jr  »  -C?  H0k..  (20) 

Equation  20  is  valid  for  an  arbitrary;  lossless  isotropic 
medium.  To  obtaia  the  analog  of  the  ray-normal 
equation  of  Landau,  ond  Lifschilz:  [1959],  we  next 
specialize  (20)  to  dielectrics.  For  this  ease  A'  «=■  w/o 
where  o  >=>  V  «  phase  velocity,  lit  .this  limit  (20) 
becomes 

dk/dt  -  k(k-Va)  -  -Vc(x„0  (2!) 

Equation  21  determines  the  motion  of  tiie  ray  normal 
in  dielectrics  which  vary  slowly  witlt  both  position 
and' time.  As  demonstrated  by  Lunduu  and  Lifschitz, 
who  obtained  Cite  same  equation  for  tiie  ease  when 
a  varies  with  position  only,  (21)  predicts  a  bendinp 
of  tiie  rays  toward  the  region  where  v  is  smaller’. 

We  also  note  from  (20)  and  (21)  that  if  o  is  hide- 
rpendent  ofposilio.tvbuUlocs  depend  on  timcf  w«  have- 

ilk/ilt  «  0  (22) 

Therefore,  ns  expected,  the  ray  does  not  change  its 
direction  of  propagation  in  media  which  vary'  only 
witlt  time, 

2.2.4.  Temporal  discontinuities.  It  is  often  desirable 
to  know  the  behavior  of  w  and  k  when  tiie  properties 
of  the  medium  are  suddenly  altered.  For  example, 
suppose  We  have  a  dielectric  In  which  <  »  e,(x)  for 
/  <  f,„niul  t  »  «j(x)  for  /  >  /,.  To  study  the  behavior 
of «  and  k  when  temporal  discontinuities  occur  let  us 
integrate  (4)  from  I,  —  8  to  /,  +  5  (note  that  (4) 
is  noUtrictly  valid  for  $  «  0).  We  obtain 

k(x,  f,  *  5)  -  A(x,  t,  ~  8)  «  -  f" (Vie)  fit  (23) 

In  tiie  limit  ns  8  0  tiie  right-hand  side  of  (23) 

vanishes,  unless  Vu>  lias  a  delta-function  behavior; 
therefore, 

k(x,  t,  -F  6)  «  k(x,  /,  -  8)  (24) 

-Since  (24)  implies  that  both  the  magnitude  and  direc¬ 
tion  of  k  cannot  change  instantaneously  then  from  (3) 
we  may  write 

A'.Mx.  1 1  T  8),  x,  1 1  I  v] 

-  A'  Mx,  c,  -  8),  x,  I,  -  8]  (25) 

For  n  dvkctrie.  in  which  K  -  ±«fx,  t%u,i(x, !)]'  \ 
(25)  yields 

utx,  t  +  8)  H<d*l'Mx)f  Vfx.  I  -  8)  (26) 

The  po-iihc  sign  in  (26)  is  oppropnatc  for  the  wave 
traveling  along  k,  while  the  negative  sign  is  appro- 
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priate  for  the  reflected  component  which  travels 
along  —  k.  This  latter  component  is  negligible  in 
the  WKB  limit.  In  the  limit  of  a  spatially  homogeneous 
dielectric,  (26)  reduces  to  the  previous  result  of 
Morgen  tliahr  [1958). 

2.3.  Velocity  of  Propagation  of  Values  of  k  and  y 

Onc  of  the  points  we  have  made  is  that  in  a  space* 
time  varying  medium  neither  y  nor  k  is  an  invariant 
as  one  moves  with  the  group  velocity  V  «=  Viu. 
However,  it  is  possible  to  define  a  new  velocity  Vu> 
such  that  k  will  be  constant  when  the  observer  moves 
With  this  velocity.  To  obtain  V'*1,  let  us  suppose  that 


the  solution  of  <5)  if  given  by  k(x,  /).  Then 
define  a  velocity  V'  through  the  equation 

wc  can 

(V'-V)k  «*  (V  W)k,t 

(27) 

IfJhis  wcreJonc,  (5)  could  be  rewritten  as 

Ok/Ot  [(V  H-  VO’Vlk  «  0 

(28) 

’from  which  wc  immediately  identify 

V<0  «  V  4-  V' 

(29) 

Equation  27  can  be  solved  for  V';  by  standard  matrix 
methods.  It  is  interesting  toconsiderthe  one  dimen¬ 
sional  |imit  (i.c„  O/Ox  ®  0/0y  »  0).  Then  (27) 
becomes 

n.v‘H§)(  «  (OU'/Oz),', 

(30) 

which  is  -readily  solved  for  V\  Using  tins  result, 
along  with  (6),  in  (29)  then  gives 

^  -  (f L + (f ).„/(D, 

(31) 

Therefore,  we  have  demonstrated  it  is  possible  to 
define  a  velocity  V'k)  such  that  the  observer  moving 
with  this  velocity  secs  constant  values  of  k.  A  similar 
Argument  holds  for  y. 

2,4;  Approximate  Solution  of  Equation  8 

In  many  instances  it  is  difficult  to  obtain  exact 
solutions  of  (5)  and  (8)  for  space-time  varying  media. 
This  is  especially  true  when  dispersion  is  present. 
However,  when  the  frequency  shift  in  propagating 
through  the  medium  is  small  compared  with  the 
transmitter  frequency,  it  is  possible  to  solve  (8)  by 
iteration,  Let  us  consider  the  limit  when 

w*1  J  (0  K/Ot)  (Is  <.<  I  (32) 

wlterc,  as  before,  k  «=>  A'(y,  x,  /)  and  ds  »  fc-c/x. 


Tlie  integral  is  along  the  ray  path  in  the  medium 
between  the  transmitter  and  the  observer.  For  the 
special  case  of  a  dielectric,  this  requires  that,  in 
addition  tp  the  condition  that  the  medium  vary  slowly, 
the  path  length  in  the  medium  cannot  be  so  large 
that  J  A'  ds  »  I.  When  (32)  is  satisfied,  we  may  neglect 
the  right-hand  side  of  (8),  and  conclude  as  a  lowest- 
order  approximation  that  y  is  a  constant  of  the 
motion.  Furthermore,  if  (as  is  usually  the  case)  we 
specify  that  «  =  «„  for  ull  time  at  the  position  of  the 
transmitter,  then  y  y0  for  all  x  and  /  (for  which  (32) 
is  still  satisfied).  Using  y  a*  co0  then  gives 

|V<5  i=»  f(wo/c )«(w0.  x,  t)  (33) 

where  n(u,  x,  t)  is  the  index  of  refraction,  defined  by 
|kj  *>  (u/c)n.  Equation  33  is,  of  course,  the  standard 
eikonal  used  in  ray  optics.  Similarly,  for  0  we  have, 
substituting  |k[  =■  (y0/c)ji(y0,,x,  f), 

L/c  J  «(«0,  x,  r)k- c/x  —  ayj  (34) 

where,  as  before,  k  is  the  unit  vector  normal  to  the 
phase  surface,  and  is  determined'  by  solving  (33). 
Finally,  since  w  -=>  — 0(j>/0t  wc  have  for  .the  first 
iteration  to  the  instantaneous  frequency 

«  cx  y„  —  (y0/c)  J  {[<)/i(<vo,  'x,  /)J/0/)  ds  (35) 

where  r/s  =  fc-r/x  «=  path  length  along  the  ray. 
Equation  35  is  the  result  used  by  ionospheric  re- 
searoucrsvin  studying  The  poppler  shift  through’  an 
ionospheric  region  wiiicit  varies  slowly  witii  path 
position  and  time  [Wcokcf,  1958;  Kelso,  I960,  1964; 
Ginzburg,  1964;  Pennell,  1967).  It  is  evident,  from  (32) 
and  (35)  that  (35)  is  valid  only  when  the  floppier 
shift  is  small  compared  with  the  transmitter  frequency. 
For  most  problems  of  propagation  through  the 
earth’s  ionosphere,  (35)  is  an  adequate  approximation 
for  the  instantaneous  frequency,  However,  there  arc 
laboratory  plasmas  and  sonic  planetary  nlmosphercs 
(c.g.,  Jupiter)  where  (35)  may  not  be  a  good  approxi¬ 
mation,  In  addition,  the  constraint  of  (32)  may  not 
hold  in  many  spacc-tinie  varying  dielectrics,  lift  the 
next  section,  wc  will  study  the  exact  solutions  of 
(5)  and  (8)  in  some  dielectric  materials. 

3.  APPLICATION  TO  ISOTROPIC, 

LOSSLESS  DIELECTRICS 

3. 1 .  Calculation  of  the  Invariants,  y,  and  k  in  Dielectrics 

Wc  will  now  use  the  results  of  section  2  to  study  the 
propagation  in  lossless  dielectrics  with  permittivity 
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varying  slowly  with  position  and  time.  To  simplify 
the  problem  we  will  also  assume  that  the  propagation 
is  in  the  same  direction  as  V «,  which  we  choose  to  be 
along  the  z  axis  in  a  rectangular  coordinate  system. 
For  this  case  (8)  becomes 

ilw/iU  4-  ft'\<ka/0z)  =  — w(<y<1/)(ln  ft)  (36) 

where  /I  —  From  the  theory  presented  in 

(11)  and  (12)  we  know  that  the  solution  of  (36) 
will  have  the  form  Cx  =  «I»(Ca)  where  C,  =  /(<■>,  z,  t) 
and  C-i  —  v(<".  z,  0  are  the  particular  integrals  of 

ill  =  <l:/fr'  =  -c/w/|w[<)(ln  0)/dt)\  (37) 

We  discuss  the  solutions  of  (37),  below,  for  several 
special  cases. 

3.1.1.  f!  separable.  When  ft  is  separable  we  may 
write  ft  =  P\(z)  PAD-  In  this  case,  upon  combin¬ 
ing  the  first  and  third  members  in  (37)  we  get 

vflAO  =  C,  (38) 

Therefore,  the  observer  moving  with  the  velocity 
V  =  p~l  finds  that  n>ftAO  is  an  invariant.  Similarly 
combining  the  first  and  second  members  of  (37) 
wc  get 

/  PA:’)d:  -  f  ilr'/ IP  At')]  =  C,  (39) 

so  that  the  general  solution  of  (36)  when  ft  is  sep¬ 
arable  is 

«fc.  0  -  [/MOr'«I>[ /  0A:’)  <k  -  f  <lt7ftADJ 

(40) 

where  the  arbitrary  function,  «I>  is  determined  by 
specifying  boundary  conditions  on  <»  along  any  curve 
in  the  z  —  t  plane.  For  example,  if  ft*  did  not  de¬ 
pend  on  time,  and  wc  specified  tit  —  t«ii  at  z  =  0  for 
all  t,  then  *l>  =  constant,  and  therefore  «.»  =  con¬ 
stant,  as  would  be  expected. 

3.1.2.  Taylor  expansion  of  ft.  The  situation  when 
ft  is  separable,  studied  above,  does  not  usually 
occur  in  practical  situations.  However,  there  arc 
often  problems  in  which  it  is  appropriate  to  expand 
ft{z.  t)  in  Taylor  scries  in  cither  z  or  t.  For  ex¬ 
ample,  we  can  consider  the  case  in  which  wc  desire 
to  study  the  propagation  only  over  the  time  interval 
t\  <  t  <  /■..  In  that  interval  wc  may  expand  ft  in 
Taylor  series  in  t  as 

(41) 


For  the  variation  of  (41)  the  constants  Cx  and 
C»  are 

w  exp  [Sf:)]  =  C,  (42) 

l  exp  [-S(r)l  =  fdz'  PA:’)  exp  [-S(r'))  -  C,  (43) 

where  S(z )  =  /*  ftt(z)  dz.  Upon  using  (42)  and 
(43)  in  the  general  solution  C\  —  'I'fCa)  wc  get 

to(r,  f)  =  exp  [—  S(:))-I> 

•{r  exp  [-$(.-)]  -  f'  ih’ PA:')  exp  (-S(i'))}  (44) 

The  wavenumber  k(z,  t)  is  related  to  by  k  =  oft. 
Wc  can  determine  the  arbitrary  function  <!>(•*•)  by 
specifying  boundary  conditions  on  «.  For  example, 
suppose  wc  specify  that  ><>  =  «*„  for  all  /,  at  z  —  0. 
(This  condition  is  appropriate  for  the  case  of  a  plane 
wave  of  frequency  «.s,  transmitted  into  a  spacc-timc 
varying  half-space).  This  requires  that,  in  (44),  the 
function  <I>  =  constant.  Wc  therefore  obtain 

co  =  w0  exp  £  —  J  PA:’)  (45) 

In  the  limiting  case  when  ft «  and  px  arc  independent 
of  z  this  result  can  be  readily  shown  to  be  identical 
with  the  previous  result  of  Morgen  thaler  [1958]. 

It  is  also  possible  to  obtain  the  constants  of  mo¬ 
tion  for  other  variations  in  ft.  For  example,  sup¬ 
pose  ft  is  a  function  of  z  —  vj.  Then  it  can  be 
shown  that  <»(z,  t)  [1  —  v„  ft(z  —  vj)]  is  an  in¬ 
variant  of  motion,  except  in  the  limit  when  r„  ap¬ 
proaches  the  phase  velocity  in  the  unmodulated 
medium.  (This  is  known  as  the  sonic  region  and 
has  been  discussed  in  detail  by  Hcssel  and  Olincr 

U%H). 

3.2.  Time  of  Transit  in  Space-Time  Varying  Dielectrics 
Since  the  phase  and  group  velocities  in  a  spacc-timc 
varying  medium  are  functions  of  position  and  time, 
it  is  not  immediately  evident  how  long  it  would  take 
for  a  disturbance  to  travel  a  distance  L  along  a  ray. 
To  study  this  problem  let  us  consider  the  motion  of 
the  point  at  which  the  phase  <f>  =  0.  In  particular,  let 
us  suppose  <t>  =  0  for  t  =  /„  at  some  point  x„  on  a 
given  ray.  Then  the  time  t,,  at  which  <t>  =  0  will  reach 
another  point  x,.  along  the  ray  is  a  solution  of 

/  (k’</x  —  w  ill)  =  0  (46) 

To  study  the  solution  of  (46)  let  us  specialize  to  the 
case  when  k  and  V«  lie  along  the  c  axis.  If  wc  assume 


p(z,  t)  -  PA:)  +  l PA:)  +  ••• 
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that  the  initial  point  x„  is  z  =  0  and  the  observation 
point  is  at  z  =  L  we  can  rewrite  (46)  as 

Ml-.  It) 

/  ( k  dz  -  w  dt)  =  0  (47) 

J  (ii.  i.> 

Since  the  line  integral  in  (47)  is  independent  of  the 
path,  it  can  be  taken  along  any  curve  joining  (0,  /„) 
to  (/„,  t,/)  in  the  z-t  plane.  For  instance,  we  shall 
find  it  convenient  to  write 

f  k(z,  0)  dz  =  f  oi(L.  t )  dt  =  0  (48) 

v  0  •'I, 

Equation  48  is  an  integral  equation  to  be  solved  for 
the  transit  time  (/,,  —  r„).  To  understand  the  mean¬ 
ing  of  (48)  let  us  first  suppose  that  0(z,  t)  is  inde¬ 
pendent  of  t,  and  that  at  «>(z  =  0)  is  specified  to 
be  equal  to  <*,  for  all  t.  Then  since  <«(z  =  0)  = 
*>(z  =  L)  =  (48)  becomes  the  obvious  result 


U.  -  t,  -  f’  dz/  V(z)  (49) 

v  0 

where  V(z)  =  f/<„  »(z)|,/-.  In  the  other  limit  when 
p(z,  t )  is  independent  of  z,  we  get  (assuming  k(t  = 
t„)  =  k„  for  all  z)  that  tt.  is  a  solution  of 


L  = 


F(f')  dt' 


(50) 


which  is  the  result  obtained  previously  [Fantc,  1971]. 

As  an  example  of  the  application  of  (47)  to 
dielectrics  which  vary  with  both  position  and  time, 
let  us  consider  the  ease  when  fi  is  given  by  (41), 
along  with  the  boundary  condition  that  <u  =  at 
z  =  0,  for  all  t.  Using  (45)  in  (48)  then  gives  for 
the  transit  time 


At  =  t,.  -  t. 

=  exp  [S(L)J  f  0 „(:')  exp  [— S(r')l  dz’  (5!) 
Jo 

where 


S(z)  =  [  0,«) 

Jo 

3.3.  Discussion  of  Transmission  Through  a  Dielec¬ 
tric  Slab 

The  results  of  sections  3. 1  and  3.2  can  be  applied 
to  consider  the  transmission  of  plane  waves  through 
a  lossless  dielectric  slab  of  thickness  L  in  which  the 
permittivity  varies  slowly  with  space  and  time.  Let 
us  denote  the  solutions  in  the  slab  by 

E  =»  e,(z,  f,  z,,  t .)  exp  f"/  J  (u  dt  ±  k  </;)!  (52) 


where  «.»  and  k  arc  the  appropriate  solutions  of  (36) 
and  £'„(:,  /)  is  the  appropriate  solution  of  (69).  For 
example,  if  /J(z,  l)  =  fh(z)  ft>(t)  then 

£*.(-,  t ,  Zu,  1  /) 

=  t'.(0)|  Ifl,  (.-„)]/[0,  (.-)]! ,/J  |  WAO1/10AO}}1"  (53) 

Now  suppose  the  slab  occupies  the  region  0  <  z  <  L, 
and  a  plane  wave  £,  =  exp  [/w„(/  —  z/c)]  is  normally 
incident  upon  the  slab  from  the  region  c  <  0.  Then 
the  field  transmitted  through  the  slab  at  time  t„  into 
the  region  z  >  L,  will  consist  of  a  number  of  com¬ 
ponents.  First,  there  will  be  a  wave  which  (from  z  <  0) 
crosses  the  z  =  0  boundary  at  time  r,  and  is  trans¬ 
mitted  directly  through  the  dielectric,  arriving  at 
z  -  V  at  the  time  t..  Next,  there  is  a  wave  which 
crosses  the  z  —  0  boundary  at  time  t3,  and  arrives 
at  z  =  L *  at  the  time  t,  after  being  internally  reflected 
at  time  by  the  z  -  L  boundary  and  at  time  r,  by 
the  z  =  0  boundary.  Next,  there  is  a  wave  which 
(from  r  <  0)  crosses  the  z  =  0  boundary  at  time  ts 
and  arrives  at  r  =  V  at  the  time  t„  after  being 
reflected  twice  at  the  z  =  0  boundary  (at  times  r, 
and  Tj)  and  twice  by  the  z  -  L  boundary  (at  times 
r,  and  r4),  etc.  Let  us  define  R( 0,  r)  as  the  internal 
reflection  coefficient  at  r  =  O'  boundary  at  time  r, 
R(L,  r)  as  the  internal  reflection  coefficient  at  r  --  LT 
at  time  r  (for  the  ease  in  which  the  medium  is  spatially 
homogeneous  11(0,  r)  *•  R(L,  t)  «  |[t(r)/«„],/'  —  1| 
!(«(r)/ +  ir*), 7"(r)  as  the  transmission  coeffi¬ 
cient,  from  z  -  0"  to  z  -  0’  at  time  r,  and  ’/'(r)  as  the 
transmission  coefficient  from  :  =  //  to  :  =  L*. 
Then,  the  transmitted  field  at  r  =  V  can  be  written  as 

E(z  ~  C  ,t  ~  t.)  -  f(t,)lT(T,)A(rt)  exp  (/o)0t,) 

+  r(ra)/l(rj)/((0,  t, )«(£,,  r,)  exp  (iu„r,) 

+  T(t3) A(ts)R(0,  T,)lia.  r,)R( 0.  r,)R(L.  r4) 

•  exp  (/wi>rs)  L  •  •  •  ]  (54) 

where  A(r)  =  <•„(;  =  V ,  t  »  z.  =  0,  t.  -  r). 
From  (54)  it  is  clear  that  the  nature  of  the  transmitted 
field  will  be  known  once  the  times  r„  tj,  rj,  •  •  •  have 
been  determined.  Extending  the  discussion  of  section 
3.2  we  see  that  these  arc  solutions  of  the  equations 


f  k(z,  t.)  dz  = 

f  ’  w(o,  r,  dt’ 

(55) 

Jo 

Jr, 

r" 

r  r< 

-  j'  k(z,  r.)  dz  = 

J  u(L,  t')  dt’ 

(56) 

M. 

M, 

/  A(r,  r,)  dz  ■= 
Jo  . 

/  «(0,  »')  dt’ 

Jr, 

(57) 

■*w«s*j**K 


•s>\‘ 


-V^sv^\Sv , 
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Therefore  r(,  tj,  and  tj  could  be  obtained  by  first 
solving  (55)  for  r,.  Then,  using  the  solution  of  (55) 
in  (56),  the  resulting  equation  could  be  solved  for  t7. 

Hquation_54  can  also  be  applied  to  the  case  of 
transmission  into  an  infinite  half-space.  We  can 
obtain  this  limit  by  soiling  f(t,)  =  1  and  h'L,  t)  =•  0. 
We  therefore  obtain  for  Has  field  at  a  point  z.  at  time  i, 

E(z„  t.)  =  '(t.MO-,)  exp  (/«0t,)  (58) 

where  r1  is  the  solution  of  the  equation 


When  fs(z,  t )  is  given  by  (41),  we  can  use  the 
results  of  section  3.2  to  express  n  as 


k(z',t.)dz'  «  f  to(0,  t')  dt'  (59) 


a.W)  exp  [-5(2')]  dz' 

(60) 


Once  r,  is  known,  T(  n)  can  be  obtained  by  apply¬ 
ing  the  usual  WKB  methods  to  the  profile  t(z)  = 
do1  mz)  +  t,  fh(z)]\ 

As  an  example  of  the  .application  of  the  more  gen¬ 
eral  result  of  (54)  let  us  consider  the  case  in  which 
the  dielectric  slab  is  spatially  homogeneous.  Then 
(54)  becomes 


£(z  =  L* ,  t  -  i .)  =  tf(a7Xr.M(r,)  exp  (/«„r.) 

+  T(T3)A(T3)R(Tt)R(Tj)  exp  (Ioi0t3) 

+  T(ts)  A(Ti)R(Tl)R(Ti)R(T3)R(Ti) 

•exp  (loioTi)  4-  •  •  •  J  (gj) 

where 


A(r)  =  |[t(r)]/(e(r,))|v<  (62) 

and  t,„  is  the  solution  of 

fm  v(t')dt'  =  ml  m  =  1,  2,  3,  •••  (63) 

To  examine  the  various  frequency  components  pres¬ 
ent  in  the  transmitted  wave  of  (61)  we  can  Taylor- 
expand  the  functions  rm(r).  That  is 

rjl)  =  r„‘  -f  (t  -  I.XOtJOi),,  4  ...  (64) 

where 

Wf =  [e(0]/[e(rj]  =  { U(r«,))/(«(r.)]} (65) 
Using  this  result  we  may  rewrite  (61)  as: 

/:(?  "•  L' ,  i)  =  £  D,  exp  (/»,(/  -  /,)]  (66) 


n>  =  7-(r„  _ ,  )tT(/.)j^ 

Nf-I) 

•exp  (/w0r„-i*)  II  R(r,) 


r  ou.)  ■ 

=  B  pv-,)T 

L°(r>i -i). 

From  (66)  we  sec  that  the  transmitted  signal  con¬ 
sists  of  comper:  ’its  at  the  instantaneous  frequencies 
°i.  tts.  n;,.  n,,  •  •  >.  This  interpretation  assume?  that 
Bt  is  a  s!  vty  varying  function  of  time  in  comparison 
with  exp  (/tv),  which  is  the  case  in  the  WKB 
approximation.  The  importance  of  the  frequency 
components  at  n2,  n3,  fit,  •••  in  comparison  to  that 
at  nt  will  depend  on  the  amplitude  of  the  reflection 
coefficient  R.  For  /(  1  only  fl|  will  be  significant, 

but  tor  li  near  unity  this  conclusion  is  clearly  not 
true. 

Therefore,  we  sec  that,  because  of  the  spatial 
boundaries,  the  transmitted  signal  has  components 
at  1I|,  fl»,  n3,  •  •  •  and  not  just  at  Ot  as  found  by 
Morgcnthaler  [  1 958],  who  did  not  account  for  bound¬ 
ary  effects.  For  relative  permittivities  near  unity 
the  components  at  tjfl,  •••  will  be  negligible 
compared  with  those  at  £),.  However,  for  large  rela¬ 
tive  permittivities  the  higher-order  frequency  com¬ 
ponents  will  be  significant. 

4.  RAY  TRACING  METHODS 

In  the  general  case  when  the  direction  of  propa¬ 
gation  is  not  along  V*  it  is  not  a  simple  matter  to 
solve  (5)  or  (8)  for  cither  space-time  varying  dielectrics 
or  plasmas  (in  fact,  for  plasmas  we  cannot  generally 
solve  (5)  or  (8)  even  when  Vw,  is  along  the  direction 
of  propagation).  In  many  cases  it  is  acceptable  to 
use  the  approximation  of  (35),  but  for  others  this 
may  not  be  possible.  In  such  cases  it  appears  most 
appropriate  to  approximate  the  temporal  behavior  of 
the  medium  in  a  stepwise  fashion.  For  example,  in  a 
dielectric  during  the  interval  0  <  t  <  is  the  permit¬ 
tivity  can  be  approximated  by  «  =  <,(x)  for  0  < 
i  <  6.  *  =  u(x)  for  /,</</„••.,«»  t.v(x)  for 
lv-1  <  /  <  /.v,  where  |(t,  -  e,_.)/e,|  «  1.  In  a 
plasma  we  would  approximate  the  electron  plasma 
frequency  w,(x,  t)  by:  w,  =  «„(x)  for  0  <  /  <  r„ 
«r  =  w„a(x)  for  /,  </</„•••  ,  w,  «  w,.v(x)  for 
r.v-i  <  t  <  t\.  It  is  assumed  that  |(uv,  -  w,,.,)/ 
uM|  •:<  1  for  j  =  1,  2,  •  •  •  N. 

lo  illustrate  the  method  we  will  consider  a  space¬ 
time  varying  plasma.  We  suppose  that  the  plasma 
occupies  the  half-space  >  0,  and  is  spatially  stratified 


wheie 
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in  the  r  direction  only.  We  then  assume  that  a  signal 
with  frequency  spectrum  sharply  peaked  about  «  ■»  w, 
enters  the  medium  ut  an  angle  0.  relative  to  the  :  axis. 
We  will  follow  the  progress  of  two  distinct  points 
P  and  P'  on  the  envelope  of  this  signal  (in  a  dielectric, 
P  and  P'  could  represent  two  separate  values  of  phase, 
the  progress  of  which  we  follow).  If  P  is  located  at 
z  0  when  t  •=  then  by  Snell’s  law  we  have  that 
the  ray  path  followed  by  P  in  tite  time  interval  (/„,  /,) 
is  determined  from 

[I  -  w|1,*(0)/wi*],/'  sin  0„ 

-  tl  -  w»i,(s)/wi,J,/*  sin  !>(:)  (68) 

The  point  P  moves  along  this  path  with  the  group 
velocity  |  V,(r)|  =  <{l  —  «,ia(r)/w,*Jl/8  so  that  during 
the  time  interval  t,  —  /„  the  distance  S,  travelled  along 
the  ray  by  this  point  is  the  solution  of 

I,  -  I.  -  [  (Is/ (69) 

J  0 

where  the  integration  is  along  tiic  ray  path  (in  iso¬ 
tropic  media  the  ray  and  group  paths  arc  identical). 
Consider,  also,  another  point  P’  on  the  signal  envelope. 
If  P'  is  at  :  =  0  at  t  =»  then  P'  travel  scs  the  path 
given  by  (68),  except  that  the  distance  5/  travelled 
along  the  ray  by  this  point  during  the  interval  t ,  —  I,' 
is  the  solution  of 


Now  suppose  that  at  /  —  /,  —  6(5  — *  0)  the  point  P 
is  located  at  x  =  x,  (i.c.,  x  =  x,  is  the  coordinate  of  5,) 
and  the  ray  path  at  this  point  makes  an  angle  0{  with 
the  r  axis.  We  also  assume  that  at  /  =  /,  —  5  the 
point  P'  is  located  at  x  =■  x /  at  which  point  the  ray 
makes  an  angle  0/  with  the  :  axis.  At  /«=>/,  the  plasma 
frequency  is  suddenly  changed  from  u„,(z)  to  u, 

By  virtue  of  (24),  0  cannot  change  instantaneously 
(since  k  cannot  change  in  citli  r  magnitude  or  direction 
at  finite  temporal  discontinuities  in  the  properties  of 
the  medium)  so  that  0(t  —  6)  <=  0(1  -f-  a).  Therefore, 
in  the  time  interval  /,</</,  the  ray  path  followed 
by  P  is  determined  from 

-{'-["sf]T’’in«->  (7,» 

while  the  path  followed  by  P'  is  given  by 


-  {l  -  [^/~]  sin  0(z)  (72) 

where  the  new  frequencies  u,  und  «./  associated 
with  P  and  P'  (recalling  tluit  since  P  and  P'  travel  with 
the  local  group  velocity,  the  frequency  associated 
with  these  points  rentains  constant  during  the  interval 
l,  <  l  <  t,)  are 

«»  -  W  ~  «,.*(*.)  +  m) 

»»  -  !«.*  -  +  «,,;J(:,/)]l/a  (74) 

Therefore,  in  a  space-time  varying  medium,  not  only 
will  different  portions  of  the  signal  acquire  different 
instantaneous  frequencies  (even  though  both  /’and  P1 
entered  the  medium  with  the  same  frequency),  but 
different  portions  of  the  signal  will  also  traverse 
different  ray  paths,  as  is  evident  from  (71)  and  (72). 

It  is  clear,  then,  that  in  a  space-lime  varying  plasma 
we  must  ray-trace  independently  for  each  point  on  the 
signal  envelope,  since  the  ray  trajectories  for  different 
portions  of  the  signal  arc  different  (except,  of  course, 
in  the  limit  when  the  direction  of  propagation  is 
along  V«,  or  Ve).  For  each  point  P  followed,  the 
above  procedure  can  be  repeated  continually  during 
other  time  intervals  (i.c.,  /„  <  t  <  <  t  <  l„  . . .) 

until  the  location  and  frequency  of  P  at  t  «=■  tK  has 
been  obtained.  Of  course,  if  in  any  regime  of  space 
(or  time)  we  reach  the  situation  where  w  is  dose  to 
the  WKB  method  is  no  longer  valid,  und  wc  must 
perform  a  more  careful  analysis  (sec,  c.g.,  Kelso  [1964] 
and  Ginzburg  [1964]). 

5.  COMMENTS  ON  ENERGY  FLOW 

In  a  space-time  varying  medium  wc  have  found  that 
values  of  w  and  k  arc  not  propagated  with  V  =  V»«; 
therefore,  wc  should  not  be  surprised  to  find  that 
energy  flux  docs  not  flow  with  this  velocity  either. 
To  consider  this  problem  let  us  first  discuss  lossless, 
isotropic  dielectrics.  Wc  can  then  demonstrate,  upon 
using  (I)  in  Maxwell's  equations  and  employing  the 
slowly  varying  assumption,  that  e„  satisfies 

k(V\fr*e„)  -|-  e„(V’k)  H-  2(k'V)e„ 

-  —  At»«!  C, .[(<■>«/<)/)  4-  2w(<>t/<V)l  -|-  2a'((7e„/r?/)| 

(75) 

where  \f>  =  In  t.  Wc  note  that  (75)  reduces  to  the 
results  of  section  3.1.3  in  Horn  aiul  fl/o//[l959]  in  the 
limit  when  t  does  not  depend  on  time.  We  now 
consider  the  limit  when  kxVe  *=■  0.  Wc  then  have 
from  (75)  that 
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(<//f/r)(Al/V.)  -  (A'/V.)(«vaodn  »*'*>  (76) 

where  o  «  [ji„t(r,  r)]',/#  and  t l/ch  «  d/M  o  0/0:. 
Equation  (76)  may  be  formally  solved  for  e„  by 
integrating  along  the  characteristic  to  give 


e.(z.  i )  * 


/,) 

(«7/c»r')(ln  »*'*)  dt' 


(77) 


where  /,  is  the  time  at  which  the  signal  present  at  (z,  /) 
was  located  at  r  =  0,  and  the  integral  in  (77)  is  along 
the  characteristic.  Using  (77),  we  obtain  for  the  (lux 

i  =  (<Ao,/a  k.r 


[*(o.  /,)]/(«m.)  k.(0,  /,)|* 


■exp 


[*/: 


(3/dr'Kln  v)  dt' 


(78) 


We  next  return  to  (8),  specialized  to  a  dielectric.  The 
formal  solution  for  «  is 


u  =  w(0,  /,)  exp  ^  (tl/dl'XIn  v)  dt’J  (79) 

where  the  integral  in  (79)  is  again  along  the  charac- 
tcristic.  Using  (79)  in  (78)  we  find  that 

//«*  =  [e(0,  t,)/ii.),n |  MO.  /,)lVt«S(0.  r,)]| 

=  constant  (80) 

since  /,  is  an  invariant  (i.c.,  (0/0t  +  o  0/0z)tx  =  0). 
Therefore  I/ws  is  an  invariant  of  motion  in  a  space¬ 
time  varying  dielectric,  so  that,  once  «(z,  i)  has  been 
determined,  the  energy  flux  /(z,  t)  follows  immediately 
from  (80).  For  the  ease  in  which  <  depends  on  z  only, 
we  have,  from  previous  considerations,  w  =  con¬ 
stant,  so  that  for  time-invariant  media  we  retrieve  from 
(80)  the  well-known  result  that  the  flux  /  is  invariant. 
In  the  limit  when  the  dielectric  depends  on  time,  but 
not  on  position,  then  k  ■=■  u(u.t) 1/3  *=  constant. 
Using  this  in  (80)  gives  the  result  that,  in  spatially 
homogeneous  dielectrics, «(/)/  ■  constant. 

Let  us  now  consider  the  case  of  dispersive  media 
and  attempt  to  discover  whether  a  result  similar  to 
(80)  can  be  found.  We  consider  a  lossless,  isotropic 
plasma  in  which  the  wavevcctor  k  lies  along  Vw,.  If 
we  write  the  vector  potential  A  -  A,  exp  (/$),  where 
V*A  «■  0  (radiation  gauge)  and  E  =  —  0k/dt  we 
obtain 

(d/dtXkl/t  A.)  -  ( d/01  +  V  0/0z){k'n  A,) 

-  {k'n  A.W/dtK In  Vm)  (81) 
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where  V  **  (cJk/u)  =  group  velocity,  and  c'lk'J  “ 
wJ  —  wpa.  In  obtaining  (81)  we  have  assumed  that 
the  plasma  is  electrically  neutral  so  that  V‘E  *=>  0. 
If  (81)  now  is  formally  integrated  along  the  charac¬ 
teristic  we  obtain 


where  /i  is  the  time  when  the  group  was  located 
at  z  =  0.  The  energy  flux  (Poynting  vector)  I  = 
(/c/wm.)  VcY  “  ( uk/n .)  |/f.|'i  is  given  by 


1/J 


*  _  k(0.  /i)  ,  j  ,n  „  V|J 


u/i. 


|d,(0.  /,)|a 


exp 


[j[  (0/c)r')(ln  F)rfr']  (83) 


We  note  that  in  a  plasma  I/w  is  not  an  invariant 
since  u  is  no  longer  given  by 


[/: 


exp  ((7/0/')(ln  F)  dt' 


However,  it  is  possible  to  obtain  an  approximate 
invariant  involving  the  energy.  Following  Stepanov 
[1968]  we  consider  a  pulse  with  time  duration  f,  at 
position  z,  so  short  that  «  can  be  considered  constant 
from  t  to  /  +  •?.  SVe  then  have  upon  integrating  (83) 
from  t  to  t  -f  ? 


l  I  dt  at  m.”  jt  dt'  *(0.  /,)  |  A.(0,  /,) |J 

•  exp  [jj  (0/0t")(in  F)  dt" J  (84) 


To  perform  the  integral  in  (84)  it  is  appropriate  to 
change  the  variable  of  integration  from  t'  to  /,.  This 
requires  calculating  (0t'/dt,).  The  time,  r,(z,  /).  at 
which  the  signal  present  at  (z,  t)  was  located  at  z  =  0 
is  an  invariant  of  motion,  and  therefore  satisfies 
(d/01  +  V  0/dz)t,  =  0.  If  we  differentiate  this 
equation  with  respect  to  time  we  find  that 


(d/dt  +  V  0/dzX0tx/dt)  =  KO/dt)  In  VX 0ti/dt)  (85) 
Upon  solving  (85)  for  (dtx/0t)  we  have 

{dtx/dt)  =  exp  [jj  (0/dt'X In  F)  dt' J  (86) 

Using  (86)  in  (84)  we  may  write 
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W/hu  ^  (a,h)-'  f  '  *'  dt,  k{0,  /,)  |  ,4.(0,  r)|5 

•'ll 

=  constant  (87) 

where  f,  is  the  duration  of  the  pulse  at  r,  =  0,  It  is 
Planck’s  constant,  and  W  s  j,‘*'  I  dt'  =  energy  in 
the  pulse.  Equation  87  states  that,  for  a  short  pulse, 
the  number  of  quanta  in  a  wave  packet  is  approxi- 
*  matcly  invariant  (since  /,  and  f,  arc  invariants). 

6.  SUMMARY 

We  have  studied  the  properties  of  the  WKB 
solutions  in  lossless,  isotropic  space-time  varying 
media.  It  was  found  that,  in  principle,  one  can  always 
obtain  constants  of  the  motion  which  lead  to  a 
complete  determination  of  the  frequency  and  wave- 
number,  once  appropriate  boundary  conditions  have, 
been  specified.  6r.ee  w  and  k  arc  known  one  can 
readily  study  the  transmission  through  space-time 
varying  media,  as  was  illustrated  in  section  3.3. 
However,  in  general,  tlie. constants  of  motion  arc  not 
always  easily  obtained.  In  such  problems  we  have 
shown  in  section  4  that  one  can  obtain  solutions  by 
modelling  the  medium  by  a  series  of  temporal  steps. 
That  is,  the  index  of  refraction  «(x,  /)  is  approximated 
by  «,(x)  in  0  <  t  <  /„  «2(x)  in  h  <  t  <  t2,  etc.,  and 
ray-tracing  techniques  arc  applied  during  each  time 
interval. 

To  keep  our  discussion  relatively  simple,  we  have 
avoided  considering  the  effect  of  absorption.  When 
dissipation  is  present  equation  4  must  still  be  valid, 

’  since  the  function  4>(x,  t )  is  required  to  be  unique  even 
if  losses  arc  present.  However,  when  absorption  is 
present  both  «  and  k  are  complex  numbers,  if  the 
medium  varies  with  both  position  and  time.  In  addi¬ 
tion,  the  real  and  imaginary  parts  of  the  k  vector  may 
be  in  different  directions.  When  absorption  is  present 
the  group  velocity  no  longer  has  any  physical  mean¬ 
ing.  In  fact  V  =  V»u>  may  be  complex  and  may  even 
have  a  magnitude  greater,  than  the  speed  of  light, 
even  in  the  limit  of  spatially-homogeneous,  time- 
invariant  media. 

When  the  absorption  is  small  it  can  he  included  by 
perturbation  methods.  For  example  in’,  time-invariant 
media  the  ray  paths  arc  determined  neglecting  absorp¬ 


tion.  The  absorption  is  then  included  by  multiplying 
the  field  by  exp  (—  J  k"-dx)  where  the  integral  is  along 

the  ray  path,  and  k"  is  the  imaginary  part  of  k. 
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